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Abstract 

We propose a new method to analyze fluctuations in the strength function phenomena in highly excited nuclei. 
Extending the method of multifractal analysis to the cases where the strength fluctuations do not obey power scaling 
laws, we introduce a new measure of fluctuation, called the local scaling dimension, which characterizes scaling 
behavior of the strength fluctuation as a function of energy bin width subdividing the strength function. We discuss 
properties of the new measure by applying it to a model system which simulates the doorway damping mechanism of 
giant resonances. It is found that the local scaling dimension characterizes well fluctuations and their energy scales 
of fine structures in the strength function associated with the damped collective motions. 



I. Introduction 



Collective or single-particle modes of excitation which have high excitation energies usually become damped motions 
due to coupling with background compound states with high level density. Strengths of exciting such modes form a 
broad peak (when plotted as a function of the excitation energy of a nucleus) which are spread over an interval of 
excitation energy. A typical example is the giant resonances which have a damping width of MeV order while the 
resonance energy is several to a few tens MeV depending on quantum numbers and nuclear masses Strength 
distribution plotted as a function of the excitation energy is called strength function. The strength function usually 
shows a smooth profile such as Lorentzian shape, but it also exhibits fluctuations called fine structures when the 
strength function is studied in high resolution experiments [||-^. 

An important origin of the damping of the giant resonances in heavy nuclei is the doorway coupling Q|. The 
collectivity of the excitation mode is described well in terms of a coherent superposition of Ip-lh excitations, as often 
prescribed by the RPA models. On the other hand, the damping of the mode is caused by (besides the particle 
escaping whose effect is small in heavy nuclei) coupling to other kinds of excitation modes such as 2p-2h excitations or 
Ip-lh plus low-lying vibrational modes, which are called doorway states. Although the doorway states couple further 
with more complicated states, e.g. many-particle many-hole states and have their own spreading width, the total 
damping width of the giant resonances is explained reasonably well by the doorway coupling It is argued that 
the doorway states influence also the fine structure of the strength distribution [H^. However, relation between the 
doorway states and the fine structures remains to be understood since it would depend on the spreading of doorway 
states into more complicated states, which is not very well known. 

If we look at the strength function with the finest energy scale dissolving individual energy eigenstates by assuming 
small particle escaping, another mechanism of the strength fluctuation is expected to emerge. Because of the chaotic 
nature of the compound states at the high excitation energy, strengths associated with individual eigenstates would 
show statistical fluctuation following the Porter-Thomas distribution . When the particle escaping width is larger 
than mean level spacing so that individual level peaks overlap, the Ericson fluctuation emerges P-|lO| . However 
it is natural to expect that these statistical behaviors based on the random matrix theories Iqjl l|-|13[| govern the 
strength fluctuation only at small energy scales. One should rather expect that the strength function associated 
with the damping of giant resonances involves different kinds of fluctuations coexisting at different energy scales. 
(Note also there exists an argument that experimental fine structures of particle decay strength distribution under 
a certain situation resembles large fluctuations in simulated spectra generated solely on the basis of the statistical 
Porter-Thomas fluctuation, known as the pandemonium Although this is a another possible origin of the fine 

structures, the property of the statistical fluctuation should depend on the level density of compound states under 
discussion) . Thus it is important to analyze the fluctuations of strength function in connection with the energy scales 
of the fluctuation. 

In this paper, we propose an new method to perform such an fluctuation analysis of the strength functions. With 
this method we characterize the strength fluctuation with focus on its scaling behavior by taking a similar approach 
widely adopted to analyze fluctuation phenomena having self-similar multifractal structure jTsj-^ . The multifractal 
fluctuations accompany a characteristic scaling property obeying the power laws, and they are found in various kinds 
of chaotic classical dynamical systems [p^-|l7| . Furthermore multifractal fluctuations in quantum wave functions are 
found recently for electrons in disordered matter [p^-pO[. The idea is also applied to the strength function fluctuations 
pl| , p2| , and an approximate power scaling is found in a shell model calculation of the giant dipole resonances [p2[ . 
In contrast to these multifractal analyses, however, we consider in this paper more general situations where the 
fluctuation does not necessarily follow the scaling laws, and we rather try to characterize the scaling behavior which 
may be dependent on the energy scales due to different coexisting fluctuation mechanisms. For this purpose, we 
introduce a new measure of fluctuation by extending the generalized fractal dimension used in the standard 

multifractal analysis. The new measure, which we call the local scaling dimension, characterizes the scaling property 
as a function of the energy scale of fluctuation. With use of simple examples, we will demonstrate in this paper how 
the local scaling dimension carries information on the dynamics of the damping process. 

For the analysis of the strength fluctuations, the autocorrelation function has been used widely as a measure of 
fluctuation in various contexts. It was applied for example to analyze the Ericson fiuctuation in compound nuclear 
reactions to derive the particle decay width , and also to the strength of giant resonances and the decay particle 
spectra [^jj^] . The autocorrelation function and its Fourier transform is also used to characterize energy level fluctu- 
ations in the random matrix theories and chaotic quantal systems [p5|-p8[. We will show that the present approach 
and the autocorrelation analysis are related to each other and that the local scaling analysis can be used to reveal the 
energy scales involved in the strength fluctuation. 

The present method treats in a single framework strength fluctuations in wide range of the energy scale, covering 
the statistical fluctuations at the finest scale around the level spacing order as well as those associated with the 
doorway states at the larger scale. In the present paper, therefore, we use the word fine structure to denote all kinds 
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of strength fluctuations including the doorway fluctuations, which are sometimes called the intermediate structures 
in the literature. 

The paper is organized as follows. In section ||, we describe the formalism of the fluctuation analysis and introduce 
the local scaling dimension. In addition, relation to the autocorrelation function is clarified. In section til, employing 
a schematic model which simulates the doorway damping mechanism of giant resonance, we discuss in detail how the 
local scaling dimension characterizes fluctuations in the model strength function. In particular, an emphasis is put on 
the relation between the profile of the local scaling dimension and the energy scales involved in the doorway damping 
process. Finally, section IV is devoted to conclusions. 



II. Scaling Analysis of Strength Function 
A. Local Scaling Dimension 

The strength function is expressed as 

=^5,5(S-S, + So), (1) 

i 

for exciting the nucleus with excitation energy E hy sl probe which excite a mode under consideration if we neglect 
the coupling to continuum states of escaping particles and experimental resolutions for simplicity. Here and Eq are 
the energy of discrete levels and the ground state energy of the whole nuclear many-body Hamiltonian, respectively, 
and Si denotes the strength of exciting the i-th energy level. Let us assume strengths are normalized as J^i Si = I. 

By smoothing the strength function with respect to the excitation energy E with a sufficiently large smoothing 
energy width, or by putting the strengths into energy bins with large bin width, a smooth profile (e.g. a Lorentzian 
shape) will show up. If small smoothing width or bin width is used, resultant distribution will exhibit fluctuation 
associated with fine structures of the original strength function. Thus the fluctuation of the strength function generally 
depends on the energy scale with which one measures. To evaluate the scale dependent fluctuation, we apply the 
method described below. It is an extension of the scaling analysis widely used to characterize multifractal structures 

Let us consider binned distribution of the strength function S{E) by dividing whole energy interval under consid- 
eration into L bins with length e. Strength contained in n-th bin is denoted by p„, 

2(En— th bin 

To characterize the fluctuation of the binned strength distribution {pn}, moments of the distribution are introduced. 
This leads to the so-called partition function Xmi^) deflned by 

L 

Xm{e)^J2Pn=L{p^)- (3) 

n=l 

For completeness of the definition, we need to introduce an average of Xm(e) with respect to the bin boundary. This 
is because there is still arbitrariness in Eq.(||) concerning how the boundaries of bins are chosen even if the bin width 
e is fixed. To overcome this uncertainty, we make an average over different choices of the bin boundaries. When the 
average is done with sufficiently large numbers of choices so that the boundaries cover densely the energy axis, the 
averaged partition function becomes independent of the boundaries. 

It is of special interest if the partition function has a scaling property. If we consider a trivial case where the 
strength is uniformly distributed with no fluctuation, i.e., Pn — const — 1/L oc e for all L and e, then an power low 
scaling Xm(e) oc e(™~i) holds for Xm(e) with respect to the bin width e. When the fluctuation is present, the partition 
function generally deviates from this limit. An extreme case of large fluctuation is the situation where the strength is 
concentrated in a single energy level. In this case the partition function becomes constant Xm(e) = 1 for any value of 
e since the binned strength is then p„ = 1 for only one bin, and Pn> = for the others. In other words the partition 
function Xm(e) scales with zero power. If the fluctuation shows a fractal structure, i.e., it has a self-similar structure 
against the change of the scale, the partition function shows a power scaling with power different from the trivial 
cases, like 

X™(6)o.e^S;-'(™-i), (4) 
or equivalently a linear scaling between logXm(e) and loge. 
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logX™(e)«(™-l)D:™loge. (5) 

The scaling coefficient D^^.'^*'*! called the generalized fractal dimension characterizes the multifractal The 
trivial case of uniform distribution corresponds to /J^^'^t'^i — implying a one dimensional uniform object in the E 
axis. The case of the single sharp peak corresponds to D^^^^'^i = since it is just a dot (zero dimensional object). 

It should be noted however that the partition function does not follow the power scaling law in general. The 
fluctuation of the binned strength may change when we change the bin width. This will happen when the strength 
function contains fluctuations having specific energy scales. For the case of the damping of giant resonance, fluctuation 
may reflect the doorway states such as 2p-2h configurations, for which their spreading width could be one of the 
energy scales involved. Also, when we look at very small energy scale comparable with average level spacing, strength 
fluctuation associated with individual energy levels may emerge. Therefore, it is more useful to analyze "how the 
partition function scales for different energy scales" than to seek the power scaling property. For this purpose, we 
introduce an extension of the generalized fractal dimension by defining the scaling coefficients as local linear coefficients 
of logXm(e) vs. loge at each energy scale e, or equivalently defining by 

An(e = — ■ (6 

m — 1 a log e 

We call this quantity the local scaling dimension since it is defined at each energy scale e and this is a function of e. 
Note also that if the local scaling dimension Dm{e) is constant over a long interval of the energy scale e, it implies 
presence of the multifractal structure of the strength fluctuation. In this case the local scaling dimension reduces to 
the standard generalized fractal dimension ZJ^actai (-jgfjj^gjj j^y jrjq^ 

In actual calculation of the local scaling dimension, we define it by means of finite difference under the change of a 
factor 2, 

TO — 1 log 2 

rather than the derivative in Eq. (|^). Using the finite difference, the calculation is very simple for all the moments. 
The local scaling dimension D2{e) defined by the derivative has an exact form (Q) expressed in terms of the individual 
energy levels Ei and associated strengths Si, as given in Appendix]^ However, exact calculation with use of Eqs. (||) 
and (|3^) is more complicated than the calculation using the finite difference (|^) except for the second moment local 
scaling dimension -D2(e)- Note also that the finite difference gives an effective smoothing. We will discuss difference 



between the two definitions in Sec. [II B 



B. Schematic Examples 

Let us illustrate how the local scaling dimension carries information on the fluctuation of strength function by using 
some schematic examples. 



1. Lorentzian 

As the first example, we consider the Lorentzian distribution 

This represents a typical peak profile associated with many kinds of damping processes. The Lorentzian distribution 
shows a smooth variation of the strength centered around the peak. Discreteness of the energy levels is neglected. 
Note that the Lorentzian distribution has one characteristic energy scale which is the FWHM of the distribution, F. 
The local scaling dimension Dm{e) associated with the Lorentzian is shown in Fig.|^ as a function of the binning energy 
scale e. When the bin width e is smaller than F, the local scaling dimension Dm{€) takes the value close to 1, while 
Dm{t) for e much larger than F. For e around F, Dm{t) decreases sharply with increasing e. This behavior arises 
because continuous and smooth aspect of the Lorentzian distribution is dominantly measured by the scaling analysis 
using small bin width e <C F, and the local scaling dimension goes to the limit of uniform distribution Z?m(e) = 1. 
With large bin width e 3> F, most of strength is concentrated in a single bin, thus the local scaling dimension reaches 
to the sharp peak limit Dm{e) = 0. With the bin width e being the same order of the width F of the Lorentzian, a 
transition between the two limits occurs. Thus the profile of Dm{e) relates to the characteristic energy scale F of the 
Lorentzian. In fact the value of e where D„i{e) decreases most steeply corresponds well to F as seen in Fig. |[ 
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2.. GOE strength fluctuation 



If we distinguish the strengths associated with individual energy levels of highly excited nucleus, the strengths 
fluctuate from level to level on top of a smooth profile of strength distribution. The level-by-level fluctuation originates 
from the chaotic nature of excited nuclei and is known to follow generic statistical rules described by the random 
matrix theory of the Gaussian orthogonal ensemble (GOE) |^,|ll| Let us characterize the GOE strength fluctuation 
by means of the scaling analysis. 

When the strength fluctuation is described by the GOE, the strengths Si for individual states fluctuate independently 
and their statistics follows the Porter-Thomas distribution [||J|] at large A'^tot limit (iVtot being the dimension of the 
random matrix). The distribution of energy levels Ei follows a semi-circle level density in average, and the position 
of the individual energy levels fluctuates locally around the average. We treat here only the strength fluctuation by 
neglecting effects both of the average level distribution and of local energy level fluctuations. Namely we assume the 
uniform distribution of the energy levels. As discussed in Appendix influence of the local energy level fluctuation 
on the local scaling dimension is small due to the spectral rigidity of the GOE energy levels. In this treatment, the 
energy levels are chosen to be equally spaced: Ei ~ id [i ~ 1, 2, A^tot), where d is the level spacing. The local 
scaling dimension is then given by 



where I = e/d (the bin width measured in the unit of the level spacing d) is used as the scaling parameter. We give 
its derivation in Appendix 

An example of the GOE strength function is shown in Fig. |^(a). The strengths of individual states fluctuate 
strongly, i.e., there exist levels which have as large strength as factor tens of the average value of strength. This is 
because the Porter-Thomas distribution is a distribution with large skewness. Because of the strength fluctuation, 
the binned strengths also fluctuate when the bin width e is not very large compared with the level spacing d. The 
fluctuation of the binned strength {pn} decreases as e increases since the binning operates as an averaging over levels 
included in bins. 

The local scaling dimension £)„(/), Eq. (||), is plotted in Fig. ||(b), which shows how fluctuation of the binned 
strength distribution changes as the energy scale (bin width) e changes. When we use the bins whose width is 
comparable with the level spacing (l ^ 10), the binned distribution shows many spiky peaks caused by large Si 
components. This leads to small values of Dm{e). The fluctuation of the binned distribution is still large for the bin 
width of order of ten level spacings {I ^ 10^), where Dm{e) deviate significantly from the uniform limit Dm{t) — 1- It 
is washed out mostly only with large bin width of order of lOOd, where Dm(e) « 1 is realized. It is also noted that the 
local scaling dimension Z?m(e) for the higher moment (larger m) takes significantly smaller values compared with the 
second moment (Z)2(e)) especially for the small bin widths I ^ 10^. Since taking the higher moment emphasizes large 
components of the binned distribution {p„}, the local scaling dimension iDm(e) for higher m reflects scaling property 
of those components. Small value of _Dm(e) for large m indicates that the large components are distributed sparsely 
along the E axis, reflecting the spiky distribution of the GOE strength. 

We evaluated also the local scaling dimension D,„(e) based on the finite difference definition, both by using the 
analytic evaluation Eq. (|3^) and by using numerical evaluation adopting 60 realizations of GOE. They give essentially 
the same curves as the those plotted in Fig. ||(b) based on Eq. (^), and the difference is almost invisible in this plotting 
scale. 

When we include the energy level fluctuation, the local scaling dimensions is affected for small bin width e compa- 
rable to the level spacing. However, the effect is not very large and becomes negligible for e > lOd, as discussed in 



As the third example, we consider the Poisson fluctuation. If we consider experimental spectra where events are 
counted in channels which corresponds to excitation energy bins, there always exists fluctuation in the spectra that 
originates from the counting statistics. Even if the strength distribution is completely uniform, counts in bins fluctuate 
statistically obeying the Poisson distribution. Relating the normalized 'strength' p„ to the counts {r„} in bins by 
Pn = fn/N {N being the total counts), the partition function Xm is expressed in terms of the m-th moment of the 
Poisson distribution. Using the average number of counts I = (r„) as the scale of the bin width, the local scaling 
dimension Dm is given by 




m — 1 



(9) 



Appendix ^. 



3. Poisson fluctuation 
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^ l + l' P + 31 + V P + 6P + 71 + 1' l^ + 10P + 25P + 15l + l (™-2.3>4,5). (10) 

The local scaling dimension takes the value significantly smaller than one for small / (/ < 10). It monotonically 
increases with increasing I and approaches to one. The behavior is similar to that of the GOE strength fluctuation, 
Fig.^(b), but the values of _D„i for the Poisson fluctuation is closer to Dm = 1 than the GOE fluctuation since the 
fluctuation associated with the Poisson distribution is smaller than the GOE fluctuation. 

Expression (|l^) indicates that effects of the counting statistics diminish for I larger than several tens for which the 
uniform limit Dm ~ 1 is almost achieved. We neglect the counting statistics in the other part of this paper. 

C. Relation to Autocorrelation Analysis 

One often uses the autocorrelation function to characterize fluctuations of the strength function. The autocorrelation 
function for the strength function, defined by 

C2{e)= j S{E)S{E + e)dE, (11) 
= Y,S^S,5{e-E,+E,), (12) 

quantifies the fluctuation correlation as a function of the displacement energy e. Dependence of C2(e) on the dis- 
placement energy e reveals the characteristic energy scale involved in the strength fluctuations. Since the local scaling 
dimension has a similar property as discussed above, one may expect relation between the autocorrelation function 
and the local scaling dimension. In fact, we can prove that C2(e) is related to the partition function X2(e) and the 
local scaling dimension I?2(e) for the second moment. 

To this end, we introduce an integral function of C2(e) by 

B2{e)= j C2(e')rfe' = I] '5,S',^?(e I). (13) 



Comparing Eq. (13) with the closed form expression Eq. (33) for the partition function, we find an equation 



(l + f^)X2(e)-i?2(e). (14) 



Thus, the partition function is expressed as 



X2{e)^- f B,{e')de', (15) 
£ Jo 

in terms of the integral of i?2(e)- 

Using this expression for Eq. (||), the local scaling dimension i?2(e) is written as 

D.{e) = (16) 
X2(e) 

which is expressed in terms of the integral B2{e) and the double integral X2(e) of the autocorrelation function C2(e). 

The close connection between the autocorrelation function C2(e) and the local scaling dimension D2{e) can be 
seen easily for the case of the Lorentzian distribution. For the Lorentzian distribution, the autocorrelation function 
C2(e) = ^Y^Te^ decreases monotonically with increasing e and its value changes most steeply around e ^ F, where F 
is the FWHM of the Lorentzian strength distribution. On the other hand, the local scaling dimension D2{f ) is given 
with use of Eq. (|T^) by 



^ ■log((6/F) 
tan-i f - log((e/F)2 + 1) 



D2{e) = , r u J. , .^ ^ (17) 
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This function also decreases monotonically from D2 = 1 at e = with increasing e. The point of its steepest slope, 
e = 0.80r, (or e = l.SF when plotted as a function of loge) corresponds well to the FWHM of the Lorentzian. This 
correspondence holds also for the local scaling dimensions Dm of higher moments as seen in the previous subsection. 

There are two figures of merit for the use of the scaling analysis in comparison with the autocorrelation analysis. 
Firstly, the local scaling dimension £),„(£) has smooth e dependence while the autocorrelation function C2(e) is a sum 
of the delta functions (if we do not include any smoothing procedure). This difference can be seen in Eq. (16) where 
the D2{e) is expressed in terms of the integral functions of C2(e). Secondly, the local scaling dimensions -Dm(e) for 
m > 2 carry information on the higher moments of strength fluctuation. 



III. Analysis of Doorway Damping Model 

In the following we shall discuss fluctuations of the strength function which is associated with damping phenomena 
of a specific state, e.g., a collective vibrational state, embedded in the background states in the highly excited nuclei. 
Keeping in mind the damping mechanism of the giant resonances, we herewith adopt a model system in which the 
damping of a collective state takes place through coupling to doorway states which consist of only a part of the 
background states. The strength function of the collective state exhibits not only spreading of the strength common 
for all damping phenomena but also fine structures reflecting the presence of the doorway states. We shall discuss in 
detail how the scaling analysis describes characteristics of this kind of strength fluctuations. 

A. Model 

The Hamiltonian of the model is given by 

H = UJc\c) (c| + Vdoorway + -^bg, (18) 

where |c) denotes the collective state, whose energy is cUc- The third term represents the Hamiltonian for the back- 
ground states, which include the doorway states as well as the other background states. The second term Vdoorway 
represents the coupling of the collective state |c) to the doorway states. 

Here we keep in mind the picture that shell model many-particle many-hole excitations form basis states {|/i)} 
of the background states and the interaction among the basis states causes configuration mixing. As a simplified 
treatment, we assume that i/bg has diagonal term representing the unperturbed energies of the basis states, and for 
the interaction among the basis states we employ a GOE Hamiltonian. Namely, 

ffbg = 5I^MlM>(Ml+^GOE, (19) 



Hgoe = + E + h.c). (20) 



The basis energy cj^ is given by an equidistant model, 

iVbs 



{-^ + ^^)d (A* = l,---,iVbg), (21) 



where d is the level spacing of the background states and A'bg represents the total number of backgrounds. Matrix 
elements v^^ and v^'^ (/x ^ i/) for the GOE Hamiltonian are independent Gaussian random variables with the zero 
mean and the variance satisfying {{v^i^)^) = 2((w"^)^). Because of the interaction, the background states have their 
own spreading width, denoted 7 hereafter, which can be estimated as, 

7^27r ^^ ^7 (22) 



The background Hamiltonian (19) is equivalent to the model adopted in ref. |30[] and the same as the Wigner ensemble 
random matrix model |3^. 

The second term in Eq. ( p^ ) representing the coupling between the collective state and doorway states is given by 

14oorway= ^' V^°°'\\c) + h.C.) . (23) 
/X— doorway 
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Here the interaction is present only for a part of the background /i states which are supposed to be the doorway states 
of the damping, and the summation in Eq.(^3|) runs only over the doorway states. The doorway ii states are selected 
in every L /i-states (i.e., ^ states with ^ = L, 2L, 3L, ...). Accordingly, the level spacing of doorway states is D — Ld. 
The coupling matrix elements V^°°'' are random variables taken from a Gaussian distribution with zero mean. 

By diagonalizing the Hamiltonian with use of the basis consisting of all the /i states and the collective state |c), 
we obtain for each realization of random variables the energy eigenstates \i), their energy Ei and associated strength 
Si — |(c|z)p of the collective state. The strength function of the collective state 

^(S) = ^5.(5(S-£;.0 (24) 

i 

is thus obtained. Since the collective-doorway couplings are uniform with respect to the energy of the background 
states, the strength function exhibits as a global profile the Lorentzian distribution. In fact, when we average the 
strength distribution over many realizations of the Hamiltonian and make smoothing with respect to the energy E, 
the averaged strength function S{E) becomes very close to the Lorentzian whose FWHM is given by the golden rule 
estimate 

r = 27r ^^ "j^ ■ (25) 

In this paper we do not discuss the global spreading profile since it is trivial in this model and it often has simple 
profile such as Lorentzian also in the experimental data. Instead we deal only with fine structures and fiuctuations of 
the strength function. It is then useful to remove the global smooth profile by normalizing the strength Si as 

S.^Af^iP^. (26) 

Here, S{E) and p{E) are obtained by averaging the calculated strength function, Eq.([2^), and the level density 
p{E) = 5{E — Ei) over realizations and by smoothing with respect to E with use of the Gaussian weighting and 
the Laguerre polynomials as adopted in the Strutinsky method p2[ . A/" is an overall n ormaliz ation factor to guarantee 



J2i Si = 1. As we discussed in connection with the GOE strength fiuctuation (Sec. II B 2 and Appendix M), effect 
of the energy level fluctuation on the fluctuation of the strength function is small except for very small energy scale 
less than about ten times of the level spacing d. Furthermore the level fluctuations within a small energy interval 
is described by the GOE model as far as 7 is not small. Thus we neglect the fiuctuation of the energy levels in the 
following analysis. In other words, we use the level order i as the scale of the excitation energy of the levels instead 
of the energy Ei itself. This is equivalent to assume that the energy levels are equidistantly distributed as Ei = id. 
In the following, we treat unfolded strength function S{E) — ^ ■ Si6{E — Ei) thus obtained. 

The number of background states is fixed to A'bg — 8191 so that the size of the Hamiltonian matrix has the 
dimension of 8192. The unperturbed energy of the collective state is placed at the energy center of the spectrum 
by putting ujc = 0, and we set the FWHM of the global strength distribution T ~ 2000(i so that the most of the 
strengths are located dominantly in the central region of the total spectra. For the interaction matrix elements in 
i?bg, or equivalently, the spreading width 7 of the background states, we use several different values 7 — GAd, 128d, 
256d, and 512d to see the dependence on 7. The spacing D between the doorway states is also varied as _D = 64(i, 
128d, and 256d, corresponding to the number of doorways 128, 64, 32, respectively. 

By generating random numbers for the Gaussian variables in the Hamiltonian, a realization of the Hamiltonian 
is obtained. The spectra and the strength distribution are calculated for each realization. We perform an average 
over many realizations of the Hamiltonian to obtain the features which is independent of specific realizations. For 
this purpose we generate repeatedly the random numbers in Vdoorway, but fix the coefficients in TJbg (Hgoe) since 
the results do not depend very much on realizations of Hqqe because of sufficiently large dimension A^bg. For the 
ensemble average, we perform averaging over 60 realizations, which is enough to obtain almost exact results for the 



GOE ffuctuations (Sec. |IIBj) 



The partition function is calculated for the bin widths of powers of 2 of the level spacing d, i.e., e = Id, I = 
1, 2, 4, 8, 16, 4096. Since the number of levels is finite, we need to take care of the edges of the strength distribution. 
The bins placed at either edge of the spectra have smaller bin width than those in the other part of the spectra, and 
this would affect the partition function when the bin width is not very small compared to the whole spectrum width. 
To avoid the edge effect, we assumed a periodic continuation of the spectra by imposing Si+Ntot = ^i, where A^tot is 
the total number of levels A^tot = A'bg + 1- For the calculation of the partition function, we performed averaging over 
16 bin boundary positions by sliding the bin boundaries by Z/16 (we average over I boundary positions for I < 16), 
which provides sufficient averaging. 
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An example of the strength function is shown in Fig. ||(a). Quahtatively the distribution looks a Lorentzian 
distribution as a gross profile. However, the distribution has large fluctuation around the Lorentzian shape. In 
addition to apparent level-to-level fluctuations, one can recognize clustering of strengths. The latter appears to reflect 
the distribution of the collective-doorway couplings as seen from the comparison with Fig. ||(b), which shows the 
strength distribution which is obtained by assuming pure doorway states, i.e., by neglecting the interaction between 
the doorway states and the other background states. Note however that it is very hard to identify individual doorway 
states in the strength distribution of Fig. ^(a) since in this example the spreading width of doorway states, 7, is larger 
than the spacing between the doorways, D. The strengths associated with individual doorways overlap with each 
other. Fig. ^(c) shows the distribution of the normalized strength (corresponding to Fig. |^(a)) for which we apply the 
scaling analysis. 

Fig. ^(a) shows the nearest neighbor level spacing distribution (NND) for D = 128d and 7 = 256d. The NND 
follows the Wigner distribution ||] perfectly, indicating that the fluctuation of the energy level at a short energy 
interval obeys the GOE random matrix limit. By using the same Hamiltonian as our background Hamiltonian i?bgj 
Guhr et al. showed that the GOE limit of the NND is reached as the spreading width of the background states 7 
becomes larger than several tens of the level spacing d. Fig. |4|(b) shows the statistical distribution of strength Si of 

individual level for D = 128(i and 7 = 256(i. Here the histogram of distribution of square root strength is plotted 
together with the Port er-T homas distribution drawn by the dashed curve. Note that the Porter-Thomas distribution 
P{S) = exp(— S'/2)/\/27rS' has a Gaussian shape in this plotting. Deviation from the Porter- Thomas distribution 
is clearly seen. The excess of large and small strengths compared to the Porter-Thomas distribution indicates that 
the strength fluctuation is larger than the GOE limit. The deviation is caused by the clustering of the strength 
distribution, which originates from the fine structures of the strength functions associated with doorway states of the 
damping of the collective state. 



B. Scaling Analysis 

1. Characteristics of Dm{^) 



Let us analyze the strength fluctuation by applying the scaling analysis which we introduced in the previous section. 

Figures]^ and ^ display the partition functions Xm{i) and the local scaling dimensions Dm{e), respectively, for the 
strength functions calculated with four different values of 7 and fixed value of Z? = 128d. From the partition function 
Xm(e), one sees that the fluctuation does not follow the power scaling law (or linear relation in the log- log plot). 
Correspondingly, the local scaling dimension -Dm(e) varies as a function of e. The local scaling dimensions in the four 
cases show common features. Namely, at small values of e, Dm{e) monotonically increases with e. As e increases 
further, Dm{e) starts to decrease, but it turns to increase again. At the largest value of e comparable to the whole 
spectrum range, Dm{() converges to the uniform limit Dm{() = 1- 

Comparing with the local scaling dimension Dm{e) for the GOE strength fluctuation indicated with dashed curve 
in the figure, one sees that the behavior of I?m(e) at small e follows the GOE limit whereas it deviates from the GOE 
limit as e increases. It is seen that the energy scale e where Dm(e) starts to deviate from the GOE limit, which we 
denote e*, has strong correlation to the value of 7. As we discussed in terms of the level spacing distribution, the 
GOE limit is realized as far as the small energy interval is concerned. However, mixing between the background states 
(including the doorway states) takes place only within a finite energy interval, that is, the spreading width 7 of the 
background states. Therefore, the fluctuations far below the energy scale 7 will obey the GOE limit, while those 
comparable to or larger than 7 may deviate from the GOE limit. In fact, one observes in Fig. || that the energy scale 
e* where Dm{i) starts to deviate from the GOE limit is related approximately to 7 as e* ^ 1/57. 

Figure || also shows clearly that the decrease of Dm(e) at the intermediate energy scale e depends strongly on the 
parameter 7. In fact, the value of e, denoted e**, where Dm(e) decreases most sharply corresponds well to the value 
of 7, indicated by the arrow in the figure, and it moves as one changes 7. Namely e** '--^ 7 is observed for all panels in 
the figure. The close connection between 7 and the decrease of £>m(e) can be explained as follows. The presence of 
the doorway states causes the fine structure or the clustering of the strength distribution. When the spreading width 
7 of the doorway states is small, the doorway states form many peaks (clusters) where width of individual peaks is 
given by 7. Using the bin width e comparable with the spreading width 7, the scaling analysis uncovers the profile of 
the peaks associated with the doorways. Since the profile will look the Lorentzian shape of width 7, the local scaling 
dimension -Dm(e) decreases sharply around e '--^ 7 as we found for the smooth Lorentzian distribution (Sec.ilBl). 
When 7 is larger than the spacing D of the doorway states, the peaks overlap with each other and one cannot trace 
the individual doorway states. However, the fiuctuations of the strength distribution still reflects the spreading width 
7 of doorway states since the fiuctuations associated with the doorway states are not completely smeared out. Thus 
correspondence between the spreading width 7 of doorway states and the point of steepest decrease in D„i{e) holds 
also for the cases of 7 > D. This situation corresponds to Figs. 0(c) and (d). 
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The increase of -Dm(e) at the largest energy scales arises from another features of the doorway states, that is, the 
fluctuation in the collective-doorway coupling. If we consider only the doorway states and neglect their coupling to 
other backgrounds, the strengths are concentrated only on the doorway states, as shown in Fig. ||(b). In this limit the 
strengths of the doorways follow a GOE strength fluctuation, since the collective-doorway coupling matrix elements 
is chosen as Gaussian random variables. It should be noted that the energy scale associated with this GOE is not the 
level spacing d of the whole spectrum, but the spacing D between the doorways, which is much larger than d. The 
dotted curve in Fig. ^ represents the GOE strength fluctuation characterized by the level spacing D of the doorway 
states. When the bin width e is taken larger than the spreading width 7 of the doorway states, the fluctuations at 
smaller energy scale are smeared out, and only the fluctuation of the collective-doorway coupling remains. In all 
cases in Fig. g, the local scaling dimensions follow this curve at the largest energy scales, indicating that the strength 
fluctuation in this domain reflects the structure of the collective-doorway coupling. 

One may also notice that the dip of £'m(e) around 6^7 becomes small and Dm{e) approaches to the GOE limit as 
one increase the value of 7. It is also noted that for small value of 7, the local scaling dimension Dm{f) for the higher 
moments m takes significantly smaller value around the intermediate e than Dm{e) for the lower moments (eg. vs. 
D2), while, for large 7, Dm{^) for all moments takes similar value. The higher moments are sensitive to components 
of the large strength p„. Therefore, small value of for higher moments implies that the strength distribution 

is very spiky. In other words the strengths are clustered and form large peaks. On the contrary, similar values of 
£'m(e) close to 1, obtained for large 7 values, implies relatively uniform distribution of strengths. The above two 
features characterize the way how the fine structures (or clustering) associated with the doorway states diminishes as 
7 increases. 

In Fig. !?[ we show dependence on the parameter D, the spacing of the doorway states, with the spreading width 7 
being kept constant. It is seen that the values of local scaling dimension Dm{e) at larger energy scale e depends on 
the spacing D of doorway states. When D is small (e.g. Fig. 0(a), where D = 64d), 15m (e) takes larger value close 
to the GOE limit than the case of larger D. This implies that when the strength of each doorway states overlaps 
significantly with each other, namely, the spacing D of doorway states is much smaller than spreading width 7, the 
fluctuation associated with the doorways tends to be smeared. On the other hand, when 7 Z), the fluctuation 
associated with the doorways stands up clearly, then D.„i{e) exhibits large deviation from the GOE limit. However, 
it is more important to note that the spacing D of the doorways does not affect very much the behavior of Z?„i(e) at 
smaller energy scale. In both cases of Fig. ^ the values e* and e** where I?,n(e) deviates from the GOE and decreases 
most sharply, respectively, have relations to 7, e* ~ 1/57 and e** ~ 7, as was observed in Fig.|^, indicating that these 
features hold in spite of the change of D values. 

In Fig. H we also plot with the dashed-dotted curve £'2(e) which is evaluated by means of the exact expression 
Eq. (34) based on the derivative definition Eq. (^. It is seen that difference between the exact evaluation and the 
evaluation Eq. (|^) based on the finite difference is almost indistinguishable. This justifies the use of the finite difference 
evaluation of the local scaling dimension. 



2. Dm (e) as a measure of doorway spreading width 

The close relation between the spreading width 7 and the profile of the local scaling dimension may be used as a 
tool to estimate the spreading width 7 of the doorway states. The number of doorway states which are involved in 
the damped collective states such as the giant resonances is not small. In such cases, it is difficult to evaluate the 
spreading width by distinguishing individual doorway states from the whole strength distribution since the strengths 
associated with each doorway states may overlap with each other. On the other hand, since the scaling analysis 
discussed in this paper does not assume separation of the peaks, it is applicable in such overlapping cases. Indeed, the 
numerical analysis presented above indicates that the characteristic behavior of Dm{t), namely the decrease around 
the energy scale comparable with 7, shows up even in cases where the spreading width 7 is larger than the spacing 
D between the doorway states. Thus one may estimate the spreading width 7 by locating the point e** where _Dm(e) 
decreases most steeply. 

It is instructive to compare with the autocorrelation analysis, which is often used to analyze the strength fluctuations. 
As discussed in the previous section, there is the relation between the autocorrelation analysis and the scaling analysis, 
in particular, between the autocorrelation function C2(e) and the local scaling dimension D2{e) for the second moment. 
Figure || shows the autocorrelation function C2(e) for D = I28d and 7 = 256d, which is also calculated for the 
normalized strength distribution (e.g. Fig. ^(c)) and averaged over 60 realizations as done for the scaling analysis. 
The autocorrelation function C2(e) exhibits the following property. It has a broad peak centered around e = which 
is built upon a flat 'background'. The broad peak contains the information of the fine structure of the strength 
fluctuations which arises from the doorway states. The width of the broad peak corresponds well to the spreading 
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width 7 of the doorway states. It is helpful to remind that the autocorrelation function used to analyze the Ericson 
fluctuation in cross sections of highly excited compound nuclei, where resonance overlaps with neighboring ones. In 
that case, the width of the autocorrelation function is related to the particle decay width of individual states. In our 
case, on the other hand, the width of the autocorrelation function is related to the spreading width of the doorway 
states. Comparing Fig. || and Fig. |8[ one sees close relationship between the local scaling dimension D2{e) and the 
autocorrelation function C2(e), both of which decrease sharply at the energy scale corresponding to 7. 

3. Analysis without ensemble average 

The results presented above are obtained by averaging over sixty different strength functions which are calculated 
for every realization of the random Hamiltonian. This is possible for the model systems. However, there will be a 
situation where one has only a single experimental data for a given nucleus so that one cannot perform ensemble 
average. When one applies the scaling analysis to a single given spectrum, the result will contain both the property 
which is common for the ensemble which the data belongs to and the behavior which is specific to the given realization 
of spectra. The latter can be regarded as fluctuation associated with different realizations. The fluctuation indicates 
an uncertainty when we infer the ensemble averaged properties from analysis of a single spectrum. To evaluate this 
uncertainty in the local scaling dimension Dm{t), we calculate the local scaling dimensions Dm(e) for each of 60 
realizations of the model Hamiltonian, and evaluate the average and the standard deviation in _Dm(e). The results 
are plotted in Fig. ^. We also plot in the same figure an example of Z)m(e) for an arbitrarily chosen realization. 

It is seen that the local scaling dimension calculated for a single realization shows e dependence which is qualitatively 
same as the ensemble averaged Z?„i(e)- The standard deviation from the average is smaller than the typical size of 
e dependence of Dm{e) for the cases shown in the figure. It is seen that the deviation increases as m increases, 
namely, as the higher moments are concerned. The standard deviation in Z)m(e) is estimated approximately as 
<y{Dm) ~ -s/Sm/V-^tot by assuming independent Porter-Thomas fluctuation for the individual strengths Si. Thus, 
the local scaling dimensions Dm{e) calculated for a single data are not affected very much by the fluctuation due to 
the realizations and keep information on the background mechanisms such as the spreading width 7 of the doorway 
states as far as a sufficient number of levels are included in the analysis and the low order moments are concerned. 

IV. Conclusions 

We proposed a new method to analyze fluctuation properties of the strength functions by means of the scaling 
analysis similar to the multifractal analysis. In order to treat the non-scaling fluctuations, we introduced the local 
scaling dimension Z),„(e) defined as a function of a scaling parameter e, which is the energy bin width subdividing 
the strength function. This quantity reduces to the generalized fractal dimension if a self-similar structure exists in 
the strength fluctuation. The local scaling dimension also has a close connection with the autocorrelation function of 
the strength function. 

Strength functions associated with damped excited modes such as the giant resonances may exhibit various fluc- 
tuations and flne structures arising from different mechanisms including the doorway structures of damping and the 
Porter-Thomas or GOE strength fluctuations associated with individual levels. Employing a model which mimics the 
doorway damping of the giant resonances, we applied the method to analyze the strength fluctuations associated with 
the doorway structures. The close relation between the profile of the local scaling dimension and the fiuctuations 
embedded in the strength function is studied in detail. We showed that the behavior of the local scaling dimension 
Dm{e) as a function of the binning energy scale e manifests clearly both the GOE fluctuation at the small energy 
scale and the properties associated with the doorway states. For example, the binning energy scale e** where the 
local scaling dimension r'm(e) decreases sharply corresponds to the spreading width 7 of the doorway states. This 
demonstrates that the local scaling dimension can be an useful measure to analyze the fluctuations and flne structures 
of the strength function associated with the damping modes in highly excited nuclei. 

It is in principle possible to apply the scaling analysis discussed in the paper to the experimental spectra and also 
to strength functions in realistic descriptions of the giant resonances in order to study the mechanism responsible for 
the strength fluctuations, such as the spreading width of the doorway states. The method itself is, on the other hand, 
quite general and not restricted to the giant resonances, and applicable to any strength fluctuation phenomenon in 
nuclei as well as in the other flelds of quantum physics. 



- 11 - 



ACKNOWLEDGMENTS 



The authors acknowledge helpful discussions with Kenichi Matsuyanagi, Toru Suzuki, and Shoujirou Mizutori. 
The numerical calculations are performed at the Yukawa Institute Computer Facility. The work is supported by the 
Grant-in- Aid for Scientific Research from the Ministry of Education, Science and Culture (No. 10640267). 



A. Explicit Expression of Partition Function and Local Scaling Dimension 

We derive here an explicit expression of Xmi^) given in terms of the strength Si and the energy Ei of individual 
levels. 

The quantity p„ defined in Eq. (^) can be written as 

p„ = (27) 

i 

by using a step function 

„ . s _ J 1 for X in n-th bin , , 



for X not in n-th bin 
Then, the partition function is also rewritten as 

Xm(e)= S^,---S,,^{Y,On{E.,,)---e,,{E,J)^i^, (29) 



where (• • ■)bin represents an average with respect to the bin boundary as explained in Sec. [I A. The quantity 
X]n^n(^ii) ' ' 'SiiiEi^) takes the value of either 1 or 0. Denoting the difference between the largest value and the 
smallest value among Ei-^ , • • • , by Di-^ ,,,i^ , the probability of the quantity being 1 is evaluated as (e — Dij ...i^)/e. 
In addition, the quantity (• • •)bin should be zero if Di-^,,,i„^ is larger than the bin width e. Thus, Xm(e) can be calculated 
as 

X™(e)- ^H---^»„.^(e-A,....,j '~^"-'"' , (30) 

il,...,im 



X <0 

As a result, the local scaling dimension Dm{e) can be explicitly written as 



^(-) ^{l : J n (31) 



In particular, for the case of m = 2, Eq. ( |30| ) and Eq. ( p2| ) lead to 

e - I A; - E. 



and 



respectively. 



X2{e) = YS,S,e{e-\E,-E,\)-^^ ^, (33) 



„ . , Y..,S.,S,e{e-\E,^E,\)\E,-E,\ 

j:.,S^.S,9ie-\E,-EMe-\E.~E,\)' ^ > 
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B. Local Scaling Dimensions for GOE Strength Fluctuations 

The local scaling dimension associated with the GOE strength fluctuation can be calculated as follows under the 
assumption that the energy levels are equally spaced Ei = id (i = 1,2,..., iVtot)- Since the strength Si of i-th level is 
written as Si = xf in terms of amplitudes Xi which obeys the GOE amplitude distribution with dimension A'^tot ||] 

Pgoe{xu ...,xn,J = 7r-^r(^)<5(l - J2 ^•?)' (35) 

i=l 

the probability distribution of {pn} defined by Eq. can be written as 

f{pi,...,PL) ^ J ... J dxi...dxN,^J{pi- ^ x'^)...5{pL- ^ a:;,^)PGOE(a;i, ■••,a;Ar^„J. (36) 

itElst bin ttE-L — th bin 

Considering the case where the bin width e — Id {I being even integer), that is, each bin contains I levels, and using 
the variable I as the scaling parameter instead of e, the ensemble average of the partition function can be calculated 

as 



,(/) = / • • • / ^Pi • ■ -dpL^PnfiPi, ■ ■ ■ ,Pl) 

•' n=l 

Ntot r(^)r(i(2m + 0) 



(37) 



I r(|)r(i(iVtot + 2m))' 

Thus, the local scaling dimension derived from Eq. ^ reads 

m+l-l m+l/2-1 

J„(V2/)^ (-log2+ log*- E log*)- (38) 

^ ) i=i i^i/2 



If we use Eq. (|6D by taking / as a continuous variable, the local scaling dimension can be written as 

m — 1 

L 

TO 



1 / 

i— 1 

As is seen from Eqs. (^) and (^9|), the local scaling dimension for the GOE strength fluctuation is independent of 
the matrix dimension. In deriving the above results, we did not perform an average with respect to the bin boundary 
for the partition function. This is not necessary since the ensemble average has the same effect for the case of GOE. 

In the above evaluation we neglected the fluctuation of the energy levels. The analytic expression which includes 
both the strength and the energy level fluctuations can be derived for the local scaling dimension of the second 
moment, D2{e). 

Using Eq. (p3|), we write the partition function as 



eX2(e) ^Jde'J de"9{e - \e' - e"\){e - \e' - e"\)A{e', e"), (40) 
where the function A{e' , e") is defined by 

A{e', e") = ^ S.SjSie' - E,)5{e" - E,). (41) 

Now, we consider the ensemble average of the function A{e\e"). Since strengths and energies can be averaged 
independently, and the average of Si and that of SiSj {i ^ j) do not depend on i and j, the result of ensemble 
average is given by 



A{e', e") = 52^(6' - 6")^ 5{e' - E,) + 5,5,^ <5(e' - E,)5{e" - E,). (42) 

i i^j 



Since the average level density is normalized ( unfolded) , 
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Y,S{e'-E,)^l, (43) 

i 

and using the Dyson's two-point cluster function P,f3|, 



J2 6ie' - E,)Sie" ^ E,) = 1 - Y{e' e"), (44) 

^^0 



the ensemble average of Eq. ([40^ can be written as 



eX2(e)/iVtot = Sfe + S^Sj{e' - / dx{e - \x\)Y{x)). (45) 



For the GOE with large A^^tot limit, we have ||,|3| 

~gl = oS~S, (a = 3), (46) 

and 

Y{e)^S{ef~J{e)D{e), (47) 
where the functions S'(e), J(e), and -D(e) are given by 

5(,) ^ (48) 

J(£)5|'±'S(,')-isgn(,). (50) 
Using (^5|)-(|50|), the local scaling dimension for the GOE strength fluctuation results in 

^2(.) = (51) 

where 

^/ N ? sin(7re) 2^., Si(7re) ,^ Si(7re) , 4 /"^'^ , , 

G(e)=e2 ^(i + _Si(^e)) + ^^(l ^) + / cosySiydy, (52) 

N / X SifTre) .Sifvre) 4 . , 4e,„., , sin^fTre), 4 /""^ , 
F(e) = e(e + a + l) + — -^(— -^ + -sin(7re)) (Si(27re) i-^) cosySiydy, (53) 

TT TT TT TT TTE TT^ 

Si, ^ '^dy. (54) 
One sees that when e is sufficiently larger than the level spacing, or practically when e ^ 10, Eq. (|5l] ) leads to 

i?2(e) =^ — (55) 
e + a ~ 1 



which is the same as Eq. (39) with m = 2. This indicates that the fluctuation of energy levels does not affect the 
local scaling dimension for e which is larger than ten times of level spacing. Note that the effect of the energy level 
fluctuation on D2{e) is not very large also for the region 1 < e < 10. At e = 4, Eq. (^Sf) gives D2 — 0.67 while inclusion 



of the level fluctuation (Eq. (51)) leads to D2 = 0.62. 



Equation ( psf ) can be used also for the Poisson fluctuation, for which the strength is constant Sf = SiSj, and there 
is no correlation among energy levels Y{x) ~ 0. Then, Eq.(p5[) leads to, up to a constant factor, 

X^=e + 1, (56) 
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and the local scaling dimension D2{e) for the second moment is given by 



D2{e) = (57) 



which is equivalent to Eq. ( p^ ) with m = 2. 
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(a)r=ioo (b)r=ioo 




FIG. 1. (a)The Lorentzian strength function, and (b) the corresponding local scaling dimension _Dm(e)- The curves 
correspond to Dm{e) for m = 2 to 5 from upper to lower. The arrow indicates the value of the width F of the Lorentzian. 




FIG. 2. (a) An example of the GOE strength function with matrix dimension of 2048. (b)The local scaling dimension Dm{e.) 
for the GOE strength function obtained after the ensemble average. 
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FIG. 3. (a)An example of the strength function obtained in the doorway damping model with parameters D = 128d and 
7 = 256d. (b)The same as (a), but coupling between the doorway states and the other background states is neglected. (c)The 
unfolded strength function S{E) corresponding to the original one plotted in (a). 
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FIG. 4. (a)The nearest neighbor level spacing distribution for the doorway damping model with parameters D = 128d and 
7 = 256d. Level spacings were unfolded. The dashed curve represents the Wigner distribution. (b)The statistical distribution 

of square root of unfolded strengths, , associated with the individual levels for D = 128d and 7 = 256fi. The dashed curve 
represents the Porter-Thomas distribution which becomes a Gaussian when plotted as a function of . Both results are 
obtained with the ensemble average over 60 realizations of the doorway couplings. 
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FIG. 5. The partition function Xm{e) for the doorway damping model with four different values of (a) 7 = 64d, (b) 128d, 
(c) 256d, and (d) 512d while D = 128d is fixed. The result is obtained with the ensemble average over 60 realizations of the 
doorway couplings. 



E 
Q 



0.8- 



0.6 



0.4 



0.2 



(a) D=128dY=64d 




E 
Q 



10 



100 
e/d 



1000 10000 



(c) D=128d Y=256d 




1000 10000 



(b) D=128d Y=128d 



(d) D=128d Y=512d 
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FIG. 6. The local scaling dimension -Dm(e) for the doorway damping model, corresponding to the partition functions plotted 
in Fig. ^ The arrows indicate the value of 7, which is chosen as (a) 7 = 64d, (b) 128d, (c) 256d, and (d) 512d. The dashed 
curve represents D2{e) for the GOE. The dotted curve represents the fluctuation associated with the doorway coupling, which 
is taken the same as a GOE fluctuation with level spacing D in the present model (see text). The dashed-dotted curve is D2{e) 
evaluated by using the derivative definition Eqs. m and (M) instead of the finite difference. 
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FIG. 8. The autocorrelation function C2(e) of the strength functions for the doorway damping model with D = 128d and 
7 = 256d, corresponding to Fig. ^c). The ensemble average is performed. The arrow indicates the value of 7. 
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FIG. 9. Comparison of the local scaling dimensions D2{e) and D5{e) between the results obtained with the ensemble average 
(solid curve) and those for a single realization of spectra (dotted curve), for D = 128d and (a) 7 = 64d, (b) 128d, (c) 256d, and 
(d) 512d. The bars indicate the standard deviation of D2{e) and -D6(e) associated with different realizations, which is evaluated 
by using the ensemble of the 60 spectra. 
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